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ABSTRACT 

Time-distance helioseismic measurements in surface- and deep-focus geometries for wave-paths that 
distinguish surface magnetic contributions from those due to deeper perturbations beneath a large 
sunspot are presented and analysed. Travel times showing an increased wave speed region extending 
down to about 18 Mm beneath the spot are detected in deep-focus geometry that largely avoids use 
of wave field within the spot. Direction (in- or out-going wave) and surface magnetic field (or focus 
depth) dependent changes in frequency dependence of travel times are shown and identified to be 
signatures of wave absorption and conversion in near surface layers rather than that of shallowness of 
sunspot induced perturbations. 

Subject headings: Sun: helioseismology — Sun: magnetic fields — Sun: oscillations — sunspots 



1. INTRODUCTION 

A problem of major importance in solar and stel- 
lar magnetohydrodynamics is the determination of sub- 
surface magnetic and t hermal constitution of sunspots 
(|Thomas &; Weissl I1992D. Applicat ion of time-distance 
helioseismology (|Duvall et"afl Il993f ) in 3-d tomographic 
imag in g of subsurface layers of suns pots (iDuvall et all 
1996 : iKosovichev fc Duvall I 119991 : Izhao et all 120011 . 



20061 ) provid ed a major milestone in the seismology 
of sunspots ([Thomas et all 1 19821 . 11988ft . Continuing 
developments in local helioseismology, notably helio- 
seismic holography accompanied by theoretical mod- 
elling, have since identified several "surface" contribu- 
tors to seismic mea sures, collecti v ely known as "surface 
magnetic effects" (iBraun I Il997t iSchunker et all 120051 : 
iLindsev & Br^nl l2005allbl : IZhao et all 12006) , that arise 
from interactions of p modes with surface magnetic field 
of a sunspot. However, such interactions have so far 
not been explicitly included in helioseismic inversions, 
for lack of a suitable model, and hence relative contri- 
butions of surface and sub-surface perturbations have so 
far not been estimated. 

Use of oscillation signals observed within sunspots 
is known to be the primary s ource of mos t surface 
effect s in sunspot seismology (Br aun I Il997t ILindsev! 
120061). as well as the associated observational errors 
( Rajaguru et al.l 120071) and syste matics in the analysis 
procedure ([Rajaguru et al.l l2006h . In this Letter, us- 
ing waves with a first skip travel distance sufficiently 
larger than the diameter of sunspot under study, and 
employing deep-focus wave-path geometries that largely 
avoid wavefield observed within the sunspot, we present 
new time-distance helioseismic diagnostics that contrasts 
near-surface perturbations with deeper ones. Signatures 
of a deep (up to about a depth of 18 Mm) increased wave 
speed region, obtained using waves in deep-focus geom- 
etry are presented in Section 3. Further support for this 
inference, from frequency dependences of surface- and 
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Fig. 1. — Ray-path diagram depicting the surface- and deep- focus 
measurement geometries: th e ray p aths correspond to model S of 
IChristen scn-Dalsgaard et al] 1119961) . Dotted lines correspond to 
the well-known center-annulus geometry, and solid and dashed lines 
correspond to annulus-annulus deep focus geometry with q=0.6 
and 1 (focus at the lower turning point), respectively 

deep-focus travel times, are presented in Section 4. 

2. DATA AND ANALYSIS METHOD 

A large sunspot (NOAA AR9057, diameter rj 32 Mm) 
that crossed the central m eridian on 28 June, 2 000 and 
observed by SOHO/MDI ([Scherrer et all [19951 ) in the 
full-disk resolution mode has been chosen for this study. 
We extracted a data cube of size 373 Mm 2 by 1024 
minute and have also obtained vector magnetograms of 
the same region observed by IVM instrument at Hawaii 
(Hannah Schunker, private communication). The data 
were pre-processed through the Stanford data pipeline 
system to remap and track the region under study. 

We choose waves having a first skip travel distance, 
A, of 50 Mm that is larger than the diameter of the 
sunspot. Apart from a phase speed filter that selects 
waves with horizontal phase speed around 54 km/s with 
a width (FWHM) of 21.25 km/s, Gaussian frequency fil- 
ters centered around 2.5, 3, 4 and 5 mHz with widths 
of 1 mHz are al so used to study freq uency dependences 
of travel times ([Braun fc Birchl l2006). We employ three 
distinct measurement schemes (Figure 1): (i) the tra- 
ditional center-annulus surface focus geometry (dotted 
lines), which, for the travel distance chosen, makes use of 



2 



Rajaguru 



a) j b) B_ (G) 




TABLE 1 
Annuli Radii, Focus Depths and B 3 
(A=50 Mm) 



ri(Mm) r 2 (Mm) z d {Mm) B 3U (G) 
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Fig. 2. — a). MDI continuum intensity (I c ) image of the 
sunspot (contours mark the umbral and penumbral boundaries, 
b) smoothed, by a 3 pixel width Gaussian, IVM line of sight mag- 
netogram rescaled to MDI resolution, c) and d) the reduced mag- 
netograms B s derived using Eqn.(l) for surface- and deep-focus 
(average of q=0.75,0.8 and 1 cases) geometries; see text for more 
details 

wave-field within sunspot either as a 'source' (out-going 
waves) or as a 'receiver' (in-going waves) but not both 
simultaneously, (ii) double skip center-annulus surface 
focus geometry, which completely avoids sunspot oscilla- 
tion signals for measurements over a region of radius 43 
Mm from the spot center, and (iii) an annulus-annulus 
deep focus geometry with focus depths, Zd, ranging from 
10.24 - 17.65 Mm (solid and dashed lines correspond to 
Zd = 16.85 and 17.65 Mm, respectively; details in Section 
2.1). 

Most studies have used the LOS magnetic field on 
the surface as a proxy for characterizing 'surface si g- 
nal' in the seismic measures (|Lindsev fc Brauril2005a Tlbl; 
iBraun fc"l3 irch 2006]) ■ We note that, in the presence of 
a deep perturbation extending from the surface, surface 
magnetic field can be a proxy for neither the surface ef- 
fects nor deeper perturbations in seismic signals individ- 
ually and exclusively. Since a seismic quantity, either a 
wave travel time or phase, measured at any single lo- 
cation (pixel) uses wave field observed at distant points 
that are spread over an annulus or pupil, a more appro- 
priate proxy to characterize the surface effect would be 
similarly averaged surface magnetograms: weighted con- 
volution of pupil functions, a{x, y), (or masks for 'center' 
and 'annulus' locations used to average the wave-field) 
with the magnetograms, B(x,y), given by, 

B s (x, y) = J a(x -x',y- y')B{x' ,y')w{r)dx'dy' (1) 

that gives each pixel a weighted average of surface mag- 
netic field over all the pixels from where the wave-field 
used in a travel time measurement has come from (Fig- 
ure 2). Weights, w(r), for magnetic field at a given 
location (x',y') is chosen to be the inverse of its hori- 
zontal distance, r, from the focus point (the measure- 
ment point, x, y), because the density of rays or wave- 
paths projected on to the surface go as inverse distance 
from the focus. We use such reduced or averaged magne- 
tograms, B s (x,y), to study any direct surface magnetic 
field induced variations in travel times. The magnetic 
field B(x,y) used here are LOS field strengths obtained 
from IVM vector magnetograms rescaled to MDI full-disk 
resolution. 



2.1. Deep-focus measurement scheme 

Two concentric annuli, centered around the point of 
measurement, are used in our typical deep-focus mea- 
surement scheme, similar to the original scheme pro- 
posed by (jDuvall et al.|[200lT ): a point from the inner 
annulus is cross-correlated with a diametrically oppo- 
site point in the outer annulus, and all such point-to- 
point correlations spanning the whole annuli (360 de- 
grees) are averaged. In this scheme, the travel distance 
A = r\ + T2, where r\ and r-i are the radii of inner and 
outer annuli, respectively; and, different focus depths 
are achieved by varying the ratio q = rx/r^ from to 
1, while keeping A fixed; q = is the traditional sur- 
face focus, and q = 1 is deep focus at the lower turn- 
ing point. This latter case me asurements were first per- 
formed bv iDuvall et alJ (|2001[ ). Here, we perform 8 dif- 
ferent focus depth measurements for the chosen A = 50 
Mm, including the surface focus; Table 1 provides fo- 
cus depths Zd, for the chosen r± and r2, computed us- 
ing ray theory for standard solar model (model S of 
iChristensen-Dalsgaard et alJ (|1996h ) and corresponding 
umbral averages B su =< B s > um bra- 

For deep focus measurements over the entire area of 
sunspot to be completley free of wave-field within the 
sunspot, both r\ and r-i should be greater than diameter 
of the sunspot, i.e. A should be greater than 64 Mm; 
but changes in travel times of waves with such large A, 
as measured in local methods, are too small to analyse 
the effects under study. Here, with A = 50 Mm, we aim 
at getting at least the measurements within umbral area 
to be free of wave-field within the sunspot, i.e. n and T2 
should at least be the sum of umbral (5.8 Mm) and full 
sunspot radii (16 Mm). From Table 1, we see that the 
last 3 (deeper most) annuli combinations roughly satisfy 
such a criterion (see also Figure 2). 

3. SEPARATING SURFACE EFFECTS: DEEP-FOCUS 
MEASUREMENTS 

In deep focus geometry, contribution of oscillation sig- 
nals from within the sunspot to a given measurement, 
as characterized by B Sl decreases as the focus depth in- 
creases; passage to B s «0 G (Table 1 and Figure 2) oc- 
curs mainly in the umbral region, and so we study the 
variation of umbral averaged travel times against, B su , 
the umbral averaged B s (umbral area is marked by the 
inner contour in Figure 2a, and it is about 5.8 Mm in 
radius). Figure 3 shows the variation of changes, with 
respect to a quiet Sun average, in umbral averages of 
one way (in- and out-going) and mean travel times, Srf u , 
8tq U and (5t^ 1 „, against B su . The error bars correspond to 
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Fig. 3. — Variation of umbral area averaged one way and mean 
travel times against B su for surface- and deep-focus geometries; 
see text for details. 

standard deviations of travel times over the umbral pix- 
els. Surface-focus geometry yields a value of B su = 60 G, 
and the corresponding travel times 6t? u and 6t% u are the 
two right most data points in Figure 3. Umbral average 
of half of double skip travel times, St^ u , which have only 
a mean signal and are completely free of wave field within 
the sunspot (cf. Section 2) corresponding to a zero B s , 
is given by the dotted horizontal line. 

Several important pieces of information about the sub- 
surface structure and interaction of p modes with the 
sunspot magnetic field are in order, in the variation of 
travel times against B su shown in Figure 3, once we un- 
derstand the variation of measurement geometry and de- 
pendent sensitivities to travel times. Firstly, because of 
geometry and the azimuthal averaging involved, the sen- 
sitivity to flows in deep focus measurements gradually 
decreases as z& increases and is zero when Zd is the lower 
turning point for waves. This happens independent of 
B s . On the other hand, passage towards B su = G in- 
volves wave-paths changing from in and out of sunspot 
to that of crossing it at some depth. Because of the 
above reasons, the asymmetric variation of one way travel 
times to a non-zero mean value of St^ u rj -10. sec at 
B su «0 G is accountable solely neither by material flow 
and its gradient in depth nor by effects localized near 
the surface due to predominant surface magnetic field 
interactions. An increased wave speed region extending 
down to, at least, a depth of about Zd — A/j/2, where 
Xh is the horizontal wavelength of waves with dominant 
power within the fr equency band used, is essential here . 
For the model S of IChristensen-Dalsgaard et al.l (|1996f ) 
the asymptotic relation (ray theory) Xh = R & c(r t ) / vr t 
yields a value of about 15 Mm for frequency v=T>.h mHz 
and lower turning point r t —18 Mm, and so the increased 
wave speed region should extend at least to within 7.5 
Mm from the focus depth around 18 Mm. This is broadly 
in agreement with early time-distance helioseismic in ver- 
sions for sound speed (|Kosovichev fc Duvall II1999I ). as 
far as the depth extent of positive sound speed changes 
is concerned. Secondly, attributing the largely B su in- 
dependent (hence, focus depth independent) values for 
5rf u (rj -10. sec.) solely to seismic signals from flow 
and sound speed requires a following surprising coinci- 
dence: (down)flow and sound speed perturbations both 
have negative depth gradients such that their signatures 
in Srf u (ingoing waves travelling against the flow) exactly 



cancel out, while those in 5t^ u (outgoing waves travelling 
with the flow) add up leading to the changes seen in Fig- 
ure 3. Frequency dependences of 5t? u , 8t^ u and St^ u , 
shown and discussed in the next Section, however, argue 
against such a scenario. Moreover, we find that half the 
double skip travel time 5t^ u (-20. sec.) is closer to 5t- u 
(-11. sec) than to 8t^„ (-35. sec), confirming the early 
findings of Braun (1997). This result, again, depends on 
wave frequency as shown in the next Section. 

In summary, results in Figure 3 show that use of wave- 
field observed within magnetic regions with significant 
B s leads to travel time measurements that cannot be ex- 
plained solely in terms of seismic signals due to flows and 
sound speed. However, our measurements in the limit of 
B su «0 G, where both surface magnetic and flow contri- 
butions go to zero, show a significant mean travel time 
signal indicating a faster wave speed region extending 
down to about 18 Mm with a resolution uncertainty of 
X h /2 « 7.5 Mm. 

4. ASYMMETRIC FREQUENCY DEPENDENCES OF ONE 
WAY TRAVEL TIMES 

Frequency dependence of mean travel times, over 
a sunspot reg i on, was shown and interpreted by 
iBraun fc BirehT p006) as an indication of perturbations 
largely confined to a region not deeper than a few Mm. 
Here, we show in Figure 4 the umbral averages St£ u , St^ u , 
St^u and the deep focus times 5r^ u (B s rj 0) against fre- 
quency; 8Tm U (B s rj 0) shown here are the ones averaged 
over the three deep most foci measurements (Table 1). 

Two important results concerning the sunspot - p mode 
interactions and deep structure of sunspot perturbations 
are contained in the frequency variations shown in Fig- 
ure 4. (1) There is a large asymmetry in the frequency 
dependences of £r/ u and 5t^ u : dominant frequency de- 
pendence is seen only in Srf u , while that in 5tq U is very 
small, for the chosen A of 50 Mm. A shallow sunspot 
perturbation cannot introduce such a direction depen- 
dent frequency variation; moreover, St^^Bs rj 0) are 
largely indepe ndent of freque n cy. S o, reasoning in the 
same lines as IBraun k, Birchl ([2006) rather suggests a 
deeper extension for wave-speed perturbation, provided 
we identify an independent physical cause for the fre- 
quency dependence of 6t? u , which we do in case (2) be- 
low. We note that for A smalle r than the diameter o f 
sunspot, as is mostly the case in IBraun fc Birchl ((2006), 
and also from our results not discussed here, both 5t£ u 
and 5tq U are equally frequency dependent. The asym- 
metric frequency dependence at larger A, then, is likely 
due to 'source' and 'receiver' locations being outside or 
inside, but not both together within the spot: whenever 
waves from outside the spot are 'received' inside, there is 
a large frequency dependence but no or little frequency 
dependence vice versa. (2) The mean times St^ u , which 
correspond to waves with both 'sources' and 'recievers' 
being outside but focussing at a first skip surface location 
within the spot, are not only closer to 6t^ u , but also have 
the same frequency dependence. This points to an inti- 
mate connection, helioseismically, between ingoing waves 
that survive the interaction with the sunspot to emerge 
on the other side and those 'observed' or 'received' within 
the sunspot: we interpret this as a helioseismic signature 
of well observed p mode absorption of sunspots. Fur- 
ther, as evident in Figure 4, the lower the frequency the 
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Fig. 4. — Frequency dependences of surface- and deep-focus 
travel times (one way and mean). Labels for different line types 
and colors identify the different measurements. For clarity, only 
error bars for mean travel times in surface focus (dashed line) mea- 
surement are drawn. 

smaller is the difference between St~^ u and St? u . This 
might indicate that mode conversion is more efficient at 
lower frequencies, for a given phase speed, than at higher 
frequencies or to a wave-number dependent sunspot - p 
mode interaction. 

5. DISCUSSIONS AND SUMMARY 

Reliable three dimensional tomographic imaging of the 
whole depth range of sub-surface layers of a sunspot is 
hard to achieve without having to use waves skipping at 
distances smaller than the size of a sunspot and hence 
without explicitly taking into account the magnetic field 
- p mode interactions. Since we lack a suitable model for 
such intera ctions and it has bee n quite clear from vari- 
ous studies (|Lindsev et al. ( 2007) and references therein) 
that significant seismic contributions arise due to surface 
magnetic field interactions, physical, instrumental and 
systematic, most sunspot local seismic inversions have 
suffered from unreliable depth descrimination of flow and 
sound speed structures. Here, we have attempted a com- 
pletely heliosiesmic diagnostic approach to checking the 
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extent of surface magnetic effects by way of combining of 
surface- and deep-focus time-distance helioseismic mea- 
surements that avoid oscillation signals observed within 
the sunspot. Usian an appropriate surface magnetic 
field proxy, B s , derived from magnetograms, we have 
been able to contrast the near-surface perturbations with 
deeper ones. 

Even though we have not proved that deep-focus mea- 
surements indeed have maximum sensitivities localized 
enough at the foci, with wave propagation calculations 
(which should form part of a future investigation), we 
believe that results presented in Sections 3 and 4 to- 
gether provide strong observational evidences for an ex- 
tended depth region of faster wave propagation beneath 
sunspots. Large asymmetry of one way travel times of 
waves focussing near the surface within a sunspot (i.e. 
when the surface magnetic field contribution is high) is 
also shown to be accompanied by asymmetric frequency 
dependence. We have identified such direction and sur- 
face magnetic field (or focus depth) dependent changes in 
frequency dependence of travel times to be helioseismic 
signatures of p mode absorption and mode conversion. 
The surface magnetic proxy B s that we have defined in 
Eqn.(l) depends only on how oscillation signals on the 
surface are used in a given measurement and so should 
prove useful in characterizing 'surface magnetic' effects 
in different methods of local helioseismology. 
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